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We study the dynamics of quantum entanglement between two Unruh-DeWitt detectors, one 
stationary (Alice), and another uniformly accelerating (Rob), with no direct interaction but coupled 
to a common quantum field in (3+l)D Minkowski space. We find that for all cases studied the initial 
O \ entanglement between the detectors disappears in a finite time ("sudden death"). After the moment 

of total disentanglement the correlations between the two detectors remain nonzero until late times. 
The relation between the disentanglement time and Rob's proper acceleration is observer dependent. 
The larger the acceleration is, the longer the disentanglement time in Alice's coordinate, but the 
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shorter in Rob's coordinate. 

PACS numbers: 03.65.Ud, 03.67.-a, 04.62.+V 



I. INTRODUCTION 



A causally disconnected (spacelike separated) pair of qubits or atoms on the same time slice can be quantum 
correlated and entangled. This evokes the notion of "nonlocality" as an innate feature of quantum entanglement, the 
precise meaning of which is a topic of sustained interest and some controversy. When one examines the quantum 
Q\ . entanglement across the event horizon of a black hole, the notion of "nonlocality" acquires an additional layer of 
meaning, pertaining not only to quantum correlations in ordinary (Minkowski) spacetime but also to some nontrivial 
, (global) spacetime structure. 

When the atoms are coupled with quantum fields, the situation becomes more interesting. Interaction with the 
quantum fields will induce decoherence of the atoms and affect the entanglement between them. It is known that 
■ the behavior of quantum entanglement is in general very different from decoherence. For the dynamics of quantum 
entanglement between two qubits different environmental settings could lead to very different results. (Compare, 
e -gi [ESQ)- Here we add in another dimension of consideration, that arising from nontrivial global structure of 
spacetime such as the existence of an event horizon, as in the spacetime of a black hole (Schwarzschild) or a uniformly 
accelerated detector (UAD). Specifically, how quantum entanglement between two detectors across the event horizon, 
one stationary and another uniformly accelerating, would evolve in time and how causality effects including that of 
retarded mutual influences would play out in these processes. This is an important ingredient for the establishment 
of relativistic quantum information theory. 

Alsing and Milburn [J] considered the quantum entanglement between two detectors (a quantum object with 
internal degrees of freedom), one is inertial (Alice) and the other is in relativistic motion (Bob, but when in uniform 
acceleration, they call it Rob) . Using the fidelity of the teleportation as a measure of entanglement, they claimed that 
the entanglement is degraded in noninertial frames due to the Unruh effect Q. In their treatment, both detectors are 
made of cavities, and the qubits are constructed by using "single particle excitations of the Minkowski vacuum states 
in each of the cavities" of a scalar field. However, as is pointed out by Schutzhold and Unruh @, the introduction of 
cavities alters the boundary conditions in the derivations of the Unruh effect @, 0, H] , and the evolution of the field 
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modes depends on the way the cavity is accelerated. Furthermore, if Rob's cavity is stationary in its local frame, then 
there is no particle creation inside, otherwise one has to take into account the dynamical Casimir effect @. 

Without the trouble associated with the cavity, Fuentes-Schuller and Mann consider free field modes in Minkowski 
space and assume the inertial Alice and the accelerated Rob are observers sensitive to different single modes. Suppose 
the quantum field is in a maximally entangled state of these two modes and the uniformly accelerated Rob is always 
in the right Rindler wedge, then from the negativity of the observed reduced density matrix (which is obtained by 
integrating out the mode in the left Rindler wedge originally sensitive to Rob at rest but undetectable when Rob 
is accelerated), Fuentes-Schuller and Mann claimed that the field state is always entangled if Rob's acceleration is 
finite, and because of the Unruh effect, the larger Rob's acceleration, the smaller the entanglement between the two 
modes. Note that the entanglement here is not between qubits or anything localized in Alice or Rob's hands; ft is 
the entanglement between the two modes spread in the whole Minkowski space and observed locally by Alice and 
Rob without considering the back-reaction of the field on the detector. Note also that both the fidelity in Q and the 
negativity in Q are independent of time. Thus there is no dynamics in these two earlier results. 

A. Main Features in this Problem 

In this paper we will consider a more realistic model with two Unruh-DeWitt (UD) detectors (atoms), which are 
pointlike objects with internal degrees of freedom described by harmonic oscillators, moving in a quantum field: Alice 
is at rest and Rob is in uniform acceleration. These two detectors are set to be entangled initially, while the initial 
state of the field is the Minkowski vacuum. 

The model has the advantage that it encompasses the main issues of interest yet simple enough because of its 
linearity to yield analytical solutions over a full parameter range. For the case of a single uniformly accelerated 
detector in a quantum field studied previously [ToL [Til ] . exact results are available in closed form. Here we will 
apply similar techniques to study how quantum entanglement evolves between two initially entangled UD detectors 
interacting through the quantum field. 

In our nonperturbative treatment the field will be evolving with the detectors as a combined closed system. The 
back-reactions from Alice's inertial detector and Rob's uniformly accelerated detector on the quantum field are 
automatically included in a self-consistent way. Moreover, our covariant formulation can account for fully relativistic 
effects. The following effects or features are included in our consideration: 

1) Unruh effect. In calculating the two-point functions of the detectors one can see that the uniformly accelerated 
detector (Rob ) w ould experience vacuum fluctuations different from those seen by the inertial detector (Alice). It 
was shown in 11] that in the Markovian regime the field state which looks like the vacuum for the inertial detector 
behaves like a thermal state for the uniformly accelerated one. This is the Unruh effect 

2) Causality and retardation of mutual influences. The initial information and the response from one detector to 
vacuum fluctuations of the field will propagate outwards and reach the other detector in finite time. But the influence 
on the other detector can in turn propagate back to the original one, which creates memory effects. These mutual 
influences should observe causality, and not propagate faster than the speed of light. They are contained in the higher 
order effects of quantum entanglement which can be calculated in our formulation. 

3) Observer-dependence. We learnt from the one-detector case that detectors with large acceleration have noticeable 
changes in time only around t = because significant time dilation will be seen by the field once |t| is large enough 

Therefore from the viewpoint of Alice and the field most of the time the detector at Rob's place looks frozen. 
On the other hand, in Rob's coordinate time Alice will take infinite time to reach the event horizon, and Alice also 
looks frozen around Rob's event horizon. These may alter the evolution of physical quantities in their dependence on 
Rob's acceleration. 

4) Fiducial time. While there is no simultaneity over space in our relativistic model, there exists a specified time 
slice (Minkowski time t = ) when we define the Hamiltonian and an initial state (as a direct product of a squeezed 
state of two detectors and the Minkowski vacuum of the field) . At every moment of time the physical reduced density 
matrix (RDM) for the two detectors will be obtained by integrating out the degrees of freedom of the field on that 
same time slice. One has to be cautious about whether the RDM depends on the time slicing or not. 

B. Key Issues of Interest 

The following key issues are addressed in this work: 

1) Disentanglement, "sudden death" and entanglement revival. Yu and Eberly [l[ discovered that, unlike the 
decoherence process, for two initially entangled qubits each placed in its own reservoir completely detached from 
the other, the disentanglement time can be finite, namely, quantum entanglement between these two qubits can see 
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a sudden death. Ficek and Tanas [2| as well as Anastopoulos, Shresta and Hu [3j studied the problem where the 
two qubits interact with a common electromagnetic field. The former authors while invoking the Born and Markov 
approximations find the appearance of dark periods and revivals. The latter authors treat the non-Markovian behavior 
without these approximations and find a different behavior at short distances. In particular, for weak coupling, they 
obtain analytic expressions for the dynamics of entanglement at a range of spatial separation between the two qubits, 
which cannot be obtained when the Born-Markov approximation is imposed. We wish to investigate in the particular 
setup of our model whether and how these distinct feature of "death" and/or "revival" manifest in the dynamics of 
entanglement. 

2) Entanglement in different coordinates. Following the well-known recipe, measures of entanglement such as 
logarithm negativity fl3T ] can be calculated in a new coordinate with a time slicing different from Minkowski times 
{e.g. Rindler time). We will study whether those measures of entanglement in a new coordinate can be interpreted 
as the degree of entanglement in Rob's clock (Rindler time). If yes, what is the difference between the entanglement 
dynamics in different coordinates? 

3) Spatial separation between two detectors. How does the entanglement vary with the spatial separation d between 
the two qubits? In (3+l)D (dimension) the mutual influences on mode functions are proportional to d^ 1 so it is quite 
small for large d even if the coupling is not ultraweak. Still, it is of interest to see whether the mutual influences 
suppress or enhance quantum entanglement, as compared to those from local vacuum fluctuations at each detector. 

These are some interesting new issues which will be expounded in our present study. 



C. Summary of Our Findings 



The results from our calculations show that the interaction between entangled UD detectors and the field does induce 
quantum disentanglement between the two detectors. We found that the disentanglement time is finite in all cases 
studied, namely, there is no residual entanglement at late times for two spatially separated detectors, one stationary 
and another uniformly accelerating, in (3+l)D Minkowski space. Around the moment of full disentanglement there 
may be some short-time revival of entanglement within a few periods of oscillations of the detectors (equal to the 
inverse of their natural frequency 0). But there is no entanglement generated at times much longer than 0(\/Vt). 

In the ultraweak-coupling limit, the leading-order behavior of quantum entanglement in Minkowski time is indepen- 
dent of Rob's proper acceleration a. When a gets sufficiently large, the disentanglement time from Alice's view would 
be longer for a larger a. From Rob's view, however, the larger a is the shorter the disentanglement time. Finally, in 
the strong-coupling regime, the strong impact of vacuum fluctuations experienced locally by each detector destroys 
their entanglement right after the coupling is switched on. 



D. Outline of this Paper 



This paper is organized as follows. In Sec. II we describe the setup of the problem, introduce the model and 
describe the measure of quantum entanglement we use. In Sec. Ill we present our calculations in the Heisenberg 
picture of the evolution of the operators and the two-point functions of the detectors. In Sec. IV we illustrate the 
results in the ultraweak-coupling limit and beyond. In Sec. V we present a discussion on a few key issues: (a) infinite 
disentanglement time in Markovian limit, (b) entanglement and correlation, (c) coordinate dependence, (d) detector- 
detector entanglement vs detector-field entanglement, e) the relation between the degree of quantum entanglement 
and the spatial separation of two detectors, and f) how generic the features illustrated by our results are. In Appendix 
A, we show the analytic form of the mode functions, while in Appendix B, the result of the case with two inertial 
detectors weakly coupled with a thermal bath is given for comparison. 



II. THE MODEL 



Consider two Unruh-DeWitt detectors moving in (3+I)-dimensional Minkowski space. The total action is given 
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FIG. 1: Alice is at rest along the world line z^(t) — (t, 1/6, 0,0) and Rob is uniformly accelerated along the trajectory 
Zg(r) = (a -1 sinh ar, a -1 cosh or, 0, 0), so the null hypersurface t — x 1 is Rob's event horizon. The initial state of the combined 
system is defined on the hypersurface t = r = 0, when the coupling between the detectors and the field is turned on. In this 
plot we choose b — 2.2 and a — 1 so b > 2a. 



+A / d 4 x<f>(x) 



dT A Q A (T A )6 4 (x» - z>X(t a )) + 



J dT B Q B {r B )S 4 (x" 



\(tb)) 



(1) 



where g^ v = diag(— 1, 1, 1, 1), Q A and Q B are the internal degrees of freedom of Alice and Rob, assumed to be two 
identical harmonic oscillators with mass mo = 1, bare natural frequency f2o, and the same local time resolution 
(so their cutoffs A and Ai in the two-point functions [To| are the same). t a and t b are proper times for Q A 
and Q B , respectively. The scalar field $ is assumed to be massless, and Ao is the coupling constant. Alice is 
at rest along the world line z^(t) = (t, 1/6, 0,0) (t a = t) and Rob is uniformly accelerated along the trajectory 



sinh ar, a 1 cosh ar, 0, 0) (t b = r) with proper acceleration a. For simplicity we consider the cases with 



b > 2a (see Fig. CD). 

Note that our model ((T]) is different from those for quantum Brownian motion (QBM) of two harmonic oscillators 
(2HO) in [3, d (and in in the rotating- wave approximation), where the cases considered are analogous to two 
Unruh-DeWitt detectors at the same spatial point. This is why there is no retarded mutual influence between 2HO's 
m . [ll|. Also here the spectrum of quantum field fluctuations felt by Alice and by Rob are different, while the 

vacuum fluctuations look the same for the 2HO's in [I3.[la.H^|. 

Suppose the coupling between the detectors and the field is turned on at t = r = 0, when the initial state of the 
combined system is a direct product of a quantum state | q A , q B ) for Alice and Rob's detectors Q A and Q B and the 
Minkowski vacuum | 0m ) for the field <&, namely, 



V>(0) 



q A ,q B ) <8> | M ) 



(2) 



Here | q A , q B ) is taken to be a squeezed Gaussian state with minimal uncertainty, represented in the Wigner function 

as 



p(Q A ,P A ,Q B ,P B ) 



1 



1 



^ (Qa + Qb? + \ (Qa - Qb? + ^ (Pa - Pb? 



j 2 {Pa + Pb? 



in which Q A and Q B can be entangled by properly choosing the parameters a and (3. 
Define the two-point correlation matrix V with elements 



(3) 



(4) 



PT 



AVA, where 



where TZ^ = {Q b (t) 7 P B (r),Q A (t), P A (t)), [i,v — 1,2,3,4. The partial transpose of V is V 
A = diag(l, 1, 1, —1). Starting with the Gaussian initial state ([3]), the reduced density matrix or Wigner function of 
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the two detectors is always Gaussian by virtue of the linearity of our model. Therefore Alice's detector Qa and Rob's 
detector Qb is entangled on time slice t if and only if [17[ 



E(i, t — a 1 sinh 1 at) = det 



< 



(5) 



where 



M 



10 
-10 
1 
0-10 



(6) 



is a symplectic matrix. We find that when E < the behavior of E is quite similar to the behavior of the negative 
eigenvalues which are connected to the logarithm negativity [HI . Thus the value of E itself is a good indicator of the 
degree of entanglement, at least in this specific model. Since it is relatively easy to obtain the analytic form of E(i) in 
the weak-coupling limit, we will calculate E rather than the logarithm negativity to determine the disentanglement 
time analytically. 

Below we calculate the two-point functions for matrix V. The uncertainty relation det[V + ifi.M/2] > can serve 
as a double-check 1171 1. 



III. EVOLUTION OF OPERATORS AND CORRELATORS 
A. Evolution of Operators 

In the Heisenberg picture [ijj HH , the operators evolve as 

it 



Qi(n) 



d 3 k 



(2tt) 3 
d 3 k 
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(7) 
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with i, j = A, B, ta =t,TB = t. qf\ 5,- , and are the (c-number) mode functions. The conjugate momenta 
are -Pa(£) = d t QA{t), Pb(t) = 9 t Qb{t), and fl(x) = d t ${x). The Heisenberg equations of motion for operators imply 



(a t 2 -v 2 )/W(x) = A 



rOQ rOO 

/ dtqf8\x-z A {t))+ drg^V(a!-z fl (r)) 
Jo Jo 



{dl +ng)g( +) (r i) k) = Ao/W(^(ri),k), 



" poo poo 

(d?-V 2 )f+\ X ,k) = A / dtq A + \t,k)5 i (x-z A (t))+ dTq B + \r,k)S\x-z B (r)) 

.Jo Jo 



(9) 
(10) 

(11) 
(12) 



which have the same appearance as the corresponding classical dynamical equations, /w) and /W look like classical 
fields generated by two pointlike sources at z A and z B . Solving the field equations pH]) and (fT2"|) , one obtains 

and /M related to gp' and g| + ' by the retarded Green's functions of the field. Inserting them into the equations of 

motion © and (jTTJ) one obtains the solutions of qV' and g. . However, the self-field induced by qf^ and q\ + ^ diverge 
at the positions of the two detectors, so one has to introduce cutoffs to handle them. Assuming Alice and Rob have 
the same frequency cutoffs in their local frame, one can do the same renormalization on frequency and obtain the 
effective equations of motion under the influence of field [l(| : 



(^+W + n?)«g ) (r) = ^ [4(T 2:1 ( .? (T))] & (4(r) - R{z B [r))) 



R{z b {t)) 



(13) 
(14) 
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(9, 2 + 2 7 a t + ^)gW(i,k) = A /<+W),k) + i ^i^gf^ (r-(^(t)),k) , (15) 

^ + 2 7 a r + ^)^ +) (r,k) = Ao/ (+) (^(r),k) + M g W^Z^lM q M ( z ° ( T ) - fl(z fl (r)),k) , (16) 

where 51 r = \/fl 2 + 7 2 is the renormalized frequency, 7 = A§/87r, and 

/ (+) (z,k) = e --*+*-x ) (17) 

R(x) = ^{ Xl -b^) 2 + p\ (18) 

X(x) ee yj (-UV + p 2 + a- 2 ) 2 + 4a- 2 f/F, (19) 

r.(«) - 4 ln 2^( X - W + ^ + ^)' ^ 

with U — t — x±, V = t + xi, uj = |k|, and p 2 = a: 2 . + x 2 . Here one can see that qa and gs are causally linked. 
The solutions of qA and qs 's satisfying the initial conditions 

/(+) (0, x; k) = e* kx , 9 t /(+) (0, x; k) = -luie^*, (21) 

9i A) (0)=9i B) (0) = l, %i A) (0)=9 T gi B) (0)=-^ r , (22) 

and /«(0,x) = 9 t /«(0,x) = <zj +) (0;k) = cV 3 <z5 +) (0; k) = ^ B) (0) = d t q<£\o) = q^ (0) = 0^(0) = 0, are listed 
in Appendix A. 

B. Two-Point functions of detectors 

When sandwiched by the initial state ([2]), the two-point functions split into 

(H^Hu) = (TZ^TZ l/ ) v + (Tl^TZ v ) a (23) 

where 

(Ti tl) n u ) v = -(o m |(^^ + ^^)|Om) 

= RC / (2^|j^ +)(tA " k) ^~ )( ^' k) ' (24) 



r< A >-r< B > 



i {n 2 /3- 2 Re (f) + rW) Rc (r£*> + ^)) + a 2 Re (r^ - ^ Re (r 
n; 2 [p 2 lm (r^ + )) Im (r[ A > + rf^) + f?oT 2 \^ (r^ - rf)) Im (r< A > - r< B ))] } , (25) 



where rjr = (q B \p < £' \qj[' \p*a > ) and p = d Tz qf \ Substituting Eqs. (|A1[) - (|A6[) , one obtains the above two-point 
functions straightforwardly. The only complication is that the integration over fc space in ([24]) may diverge, so one 
has to introduce additional frequency cutoffs corresponding to the time-resolution of the detector (Ai) and the time 
scale of switching on the interaction (A ) [13, El- After doing this, one has, for example, 

( QUr) > v = < Ql{r) )f +e(r-n){ Q 2 B (r) >« , (26) 

where t\ = a~ 1 \n{b/a) (see Fig. [T]), ( Q 2 B y®' is the two-point function of a single uniformly accelerated detector 

(expressions of ( Q% )|°\ ( P| )^ 0) , ( Q 2 A ) { ° } and ( Q% ) { ° ] have been listed in Appendix A of Ref. [U), while the 
higher order correction reads 

= § / r d/ e " r( r\ Sinn(r r" T ° < Q B {T),Q A (b- 1 - ^e-') )f 

47 2 T , ,e" r(r_T,) sin fi(r-r') r , ,, e -^-T") sin q/ _ r //) 

+ t^t / dr' i —i / dr" = x 

Q 2 J T1 i cosh ar> - \ J Tl I cosh ar" - \ 

( Qa^ 1 - o- 1 e- or '), Qa^ 1 - a"^-^") )f . (27) 
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FIG. 2: (Left) The behavior of E (solid curve) defined in (0 can be approximated by (|31 1) in the ultraweak-coupling limit with 
a sufficiently large. Two detectors are separable when E > (shaded zone), otherwise entangled. Here 7 = 10~ J , a — 0.1, 
b — 0.201, q = 1.1, f3 — 4.5, Ao = Ai = 100, fi = 2.3 and h = 1. By comparison the dotted curve is E with all ( .. ) set to be 
zero. Both curves will become positive at late times. (Right) Looking more closely, one finds that two detectors are separable, 
or totally disentangled, after yt » 4.44. 



Calculating the cross correlations ( 72.^(4), 72 b (t) } v (72 = P, Q) is also straightforward, though one has to be careful 
about the contours of integration of each term on the complex plane of k (see [l(| )■ Note that Ai is present only in 

( P A ) v and ( P B )v 

IV. DISENTANGLEMENT DYNAMICS 

A. ultraweak-coupling limit 

In the ultraweak-coupling limit (7A1 <C a, Q), the corrections due to the retarded mutual influences in mode 
functions (|A1|) - (|A6[) are 0(7) and suppressed, while the cross correlations ( TZa^b ) v (R- = P, Q) accounting for the 
response to the vacuum fluctuations are negligible. The two-point functions ( .. ) behave like (|B2[) - (|B8I) . and 

(Qi(*)) v « l^ 1 - 6 " 27 *)' (28) 

(Ql(r)) v - ^coth^(l-e- 2 ^), (29) 

and ( Pj(rj) ) v « ft 2 ( Q 2 (rj) ) v , j = A,B. It is straightforward to calculate S and determine the separability of 
Alice and Rob by inserting these two-point functions into ([5|). 



1. Evolution of entanglement in Alice's proper time (Minkowski time) 



At the initial moment t = 0, the initial Gaussian state ([3]) has 

Thus for all a 2 /? 2 ^ ft 2 , two detectors are entangled. Below we will see that quantum entanglement will, however, 
vanish at a finite disentanglement time ("sudden death"), which is usually different from the decoherence time scale 

7" 1 . 

For proper acceleration a sufficiently large such that r = a 1 sinh at <C t for t large enough, one has e 7 '* +r - ) w e 7 * 
and 

s ~-if!^ 2 -" 2 ' 32 > V21 ' <=»> 
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FIG. 3: When a — > 0, the behavior of E in the weak-coupling limit can be described by (|36[) . Here (a,/3) = (1.1,4.5) for the 
left plot (Za < 0) and (1.5, 0.2) for the right plot (Z4 > 0). The solid curves represent E in its totality, while the dotted curves 
represent E with all ( .. ) set to be zero. 



An example is illustrated in Fig. O The disentanglement time of this case looks infinite. However, the correction of 
the next order will shift the curve of X upward and make the disentanglement time finite. When Ai is large, such a 
correction is dominated by Ai terms in ( P 2 ) v : 

(Pi) v « ^(i-e-^ + ifryAi, (32) 

( P% > v « ^ c °th — (1 - e- 2 ^) + -fryAj (33) 
2 a 7r 

[see Eqs.(A4) and (A10) in Ref.[ll|]]. For tdE ^S> 1, this correction yields 

s - - w (B2 - q2 " S) V2 * + i£w < b2 - « v > 2 • < M) 



around t « which gives 



1 , 7rO 

tdE ~ — m — — . (35) 



2-f 7A 



In FiglU one nas 7*d£ ~ 4.44. 

Note that E is insensitive to a here. This implies that, to leading order in the ultraweak-coupling approximation, 
Unruh effect does little to the disentanglement between Alice and Rob from the view of Alice. 



2. Two inertial detectors in Minkowski vacuum 



For a — > 0, coth(7rf2/a) — ► 1, r — ► t, and Rob is also inertial. Suppose Alice and Rob are separated far enough, so 
the mutual influences can be safely ignored again. Then, in the weak-coupling limit with f2 3> 7A1 a — > 0, one has 

- [Z 8 (e-^-2e-^)+Z 4 }, (36) 



where 

z 8 = (h-a 2 n) 2 (p 2 - nn) 2 , (37) 

Z 4 = h 2 (/3 4 + a 4 n 4 + 6a 2 f3 2 n 2 ) - 2hft (0 1 + a 2 n 2 ) (h 2 + a 2 (3 2 ) . (38) 
It is clear that Z 8 > and Z 8 — Z± > 0. When Z4 > 0, the disentanglement time is clearly finite: 



2 7 V V Z, 



t dE « -i/l - I 1 ) • (39) 



J 8 
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FIG. 4: The solid curves represent E in its totality in Rob's point of view. "Sudden death" of entanglement occurs at 
"fTdE ~ 3.61, as can be seen in the enlargement on the right plot. The dotted curve represents E with all ( .. } v set to be zero. 
All parameters are the same as those in Fig. [5]except 7 = 10 -6 , a = 2 and b — 4.1. 



Indeed, one has jt dE rss 0.125 in the right plot of Fig. [3] 

When Z4 < 0, the disentanglement time looks infinite. But again, the corrections of the next order yields 



around t « tdE 3> 1 for large Ai, with 



z 2 = a 2 (f3 2 -nn) 2 + (3 2 (a 2 n-h) 2 >0. (41) 



This gives a finite disentanglement time 

2 7 Z 2 + ^Z 2 -Aa 2 (i 2 Zi ^ ' 

and jt dE 3.96 in Fig. [3] (left). 

The dotted curves in Figs. (left) and[3](left and right) are those Es with vacuum fluctuations of the field switched 
off, namely, with ( IZ^TZv ) v set to be zero. In Fig. [2] (left) it seems that vacuum fluctuations would reduce |E| 
thus suppressing the entanglement. This is not true. One can verify that, for sufficiently large t, the dashed curve 
in Figs. [H (left) will overtake the solid curve and then become positive, just like what Fig. (right) suggests. In 
the weak-coupling limit, the vacuum fluctuations of the field that the detectors see locally do not always suppress (or 
enhance) quantum entanglement beyond the disentanglement due to the dissipation of initial quantum fluctuations 
of the detectors (corresponding to the exponential decay of ( .. ) in time.) 



3. Evolution of entanglement in Rob's proper time (Rindler time) 

The hypersurface with constant Rindler time r extends to t < region in the left half of the Minkowski space (see 
Fig©- Suppose, even before the initial time slice with t = 0, the field state has been the Minkowski vacuum. Then 
the field state in Rindler time slicing is also a Gaussian state and the quantity E(f , r), with a similar definition to ([5]) 
but we let t — b^ 1 tanh clt here, can still serve as an measure of the entanglement between Alice's and Rob's detector 
in Rob's frame. 

For a sufficiently large, yt — * 7/6 <C 1 for large r. Then 



E(i = b 1 tanh ar, r) 



64a 2 /? 2 

and the disentanglement time is approximately 



(H 2 - a 2 f3 2 ) 



(1-e- 2 ^) coth 2 — - (1 



TdE 



-2 7 T\ 



(43) 



(44) 
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FIG. 5: E for Rob in high acceleration (solid curve) with a = 10 and b — 21. (a = .01, b = .021 for the dashed curve for 
comparison.) The large time dilation makes Rob appear almost frozen in the view of Alice at large t, and the disentanglement 
time is longer for larger a. Here 7 = 0.01, Ao = Ai = 50, and other parameters are the same as those in Fig. [2] 



in Rob's point of view (see Fig. [4]). This result may be interpreted as presenting a dynamical version of the statement 
"quantum entanglement is degraded in noninertial frames" [HQ. Note that in Rob's frame Alice will never cross the 
event horizon. So the quantum entanglement here is not the one "across the event horizon" . 

When a gets smaller, £ approaches those described in the previous subsections. When a — > 0, the value of the finite 
disentanglement time recovers the disentanglement time in the case of two inertial detectors. 



B. Beyond the ultraweak-coupling limit 



1. High acceleration regime 



When the proper acceleration of Rob is very large, Rob will reach a very high speed in a very short time, when the 
time dilation makes Rob appear almost frozen in the view of Alice, and the disentanglement process becomes slower 
than those in the case with Rob's acceleration smaller. In other words, in our setup, the larger acceleration Rob 
has, the longer the disentanglement time in Alice's clock (see Fig. [5]), while in Rob's frame t^e is shorter from (|44[) . 
The statement "a state which is maximally entangled in an inertial frame becomes less entangled if the observers are 
relatively accelerated" [|| is too simplistic and could be misleading here. 



2. strong- coupling regime 



When the coupling becomes larger, there will be oscillation emerging on top of the smooth curve of £ in the 
ultraweak-coupling limit. So around t^E the entanglement between the detectors of Alice and Rob may disappear 
and revive for several times, before they finally become separable forever. The duration of such kind of entanglement 
revival will not exceed the order of time scale I/O (see the upper plot of Fig. [5]) 

When the coupling gets even larger, the vacuum fluctuations will exert a strong impact on quantum entanglement. 
Usually this shortens the disentanglement time. Indeed, in Fig. [5] we find that, for 7A1 sufficiently large, the initial 
quantum entanglement is annulled right after the coupling is switched on. 

The higher order corrections to £ from retarded mutual influences of two detectors, while their amplitudes decay in 
time, are not always positive or negative. Nevertheless, the memory effect induced by mutual influences contributes 
little in the cases considered: Although the total corrections to £ from mutual influences become more obvious for a 
larger coupling, they are of 0(j 2 ) and remain small as seen in Fig. [6l Results with stronger mutual influences will 
be reported in a future work. 
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FIG. 6: E in the strong-coupling regime. The entanglement between Alice and Rob in their initial state is annulled by the 
strong impact of vacuum fluctuations (E > 0) right after the coupling is switched on. Here 7 = 0.1, Ao = Ai = 50, a = 1, 
b — 2.01, and other parameters are the same as those in Fig. [2] The dotted curve at the bottom are the E with all { .. ) v 
turned off. Numerically we find that the higher order corrections to E from mutual influences are less than 2% ~ 0(j 2 ). 



V. DISCUSSION 



A. Infinite disentanglement time in Markovian limit 

A criterion in the Markovian limit on finite disentanglement time has been offered by Yu and Eberly [l[. In some 
parameter range (a S (0, 1/3) there) the concurrence decays exponentially in time so the disentanglement time looks 
infinite. Here we have similar situations in the ultraweak-coupling limit as discussed in Sees. II V A II and II V A 2l In 
particular, the parameter Z4 in (|36|) plays an analogous role to a in [l[. While the entanglement gets a clear sudden 
death when Z4 > 0, the disentanglement time looks infinite otherwise. However, in Sees. IIV A ll and llV A 2l we showed 
that the correction from the next order to the ultraweak-coupling approximation will render the disentanglement time 
of the latter case finite. According to our results, it will be interesting to see whether the infinite disentanglement time 
in the Markovian cases that Yu and Eberly considered would become finite if one considers the corrections beyond 
the Markovian approximation. 



B. Entanglement and correlation 



When cross correlations vanish, two detectors are uncorrelated and the two-point correlation matrix V defined in 
(fJJ is block diagonalized so the Wigner function, or the reduced density matrix, for the detectors can be factorized 
into a tensor product of two Wigner functions or reduced density matrices for each detector. Now two detectors are 
simply separable, with a stronger condition than those for disentangled or separable states. 

In our model, the cross correlations ( TZa,T^-b ) always vanish as jt — > 00, while the two-point functions of each 



detector remain finite. We found | ( IZa^b ) \ 
residual entanglement at late times because 



e with e oc 7 at sufficiently large t. This implies that there is no 



\~jt — >oo 



Q 



\(0) 



1 A 



\(0) 



v(0) 



p- 

r r. 



AO) 



4 



(45) 



ft- 



which is a product of the uncertainty relations for each detector in steady state and is positive-definite if the coupling 
Ao is nonvanishing. (Note that ( Qa, Pa ) ~ e~ 7 * and ( Qb, Pb ) ~ cxp[— (a + 7)a~ 1 sinh -1 at] vanish at late times; 
for explicit expressions of other late-time two-point functions, see Appendix A in Ref. [ill].) In Minkowski time, two 
detectors become separable after a finite disentanglement time tdE, but not simply separable or uncorrelated until 
yt — > 00. 



C. time slicing and coordinate dependence 



Consider two events, one in Alice's world line at her proper time t = t, the other in Rob's world line at his proper 
time r = f , defined on the same time slice in some coordinate. Recall that the physical RDM for the two detectors 
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are obtained by integrating out the degrees of freedom of the field defined on the same hypersurface associated with 
a certain time slicing. Thus the RDM of the two detectors at the two events z^(t) and Zg(f) should always be 
associated with a specification of a time slicing scheme, so should the criterion of separability derived from the RDM. 
But there are infinitely many choices of time slices intersecting Alice and Rob's world line right at these two events. 
Does quantum entanglement of the detectors at z^(t) and Zg(f) depend on different choices of time slice? 

In the cases we considered in this paper, the answer is no. In our UD detector theory, the coordinate transformation 
from Minkowski coordinate to a new coordinate consistent with some other time slicing scheme such as Rindler time 
will not change the quadratic form of the action <[T|) , so the combined system in the new coordinate is still linear. If 
the initial time slice in the new time slicing scheme is identical to our t = hypersurface in Minkowski coordinate, 
and the combined system is starting with a Gaussian initial state such as $Z§, the quantum state of the combined 
system will remain Gaussian during the evolution in the new time slicing scheme. Hence the RDM of the detectors 
is still Gaussian and the criterion ([5]) can be applied in this new time slicing. The Gaussianity implies that both the 
RDM of the detectors and the criterion ([5]) are fully described by the two-point functions of the detectors ( TZ^, 1Z U ), 
which are independent of the choice of time slice connecting these two events. Thus the RDM and the separability 
of the Gaussian state at those two events are independent of time slicing. We believe that this is a general feature of 
quantum entanglement in relativistic quantum information theory. 

Nevertheless, in Sees II V ATI and HV A 3 1 we learnt that the entanglement dynamics and the disentanglement (proper) 
time are coordinate-dependent, because in general two events simultaneous in one (e.g. Minkowski) coordinate are 
not simultaneous in another (e.g. Rindler) coordinate, while quantum entanglement is a property of a quantum state 
of two events at the same time slice. [20( 



D. Detector-detector entanglement vs detector-field entanglement 

Quantum coherence in a detector is related to its purity, which also indicates the degree of quantum entanglement 
between that detector and the rest of the combined system (the field and the other detector) [18| ■ Since the combined 
system is a Gaussian state, the reduced density matrix for each detector is also Gaussian. So the purities of Alice and 
Rob's detectors are simply [ll| 

V, = m (46) 

P!)(Q])-(Pi,Q,) 

where j — A, B. Obviously the information contained in the cross correlations ( IZa, T^b ) is ignored in every combi- 
nation of Vj . We illustrate some examples in Fig. \7\ One can see that the evolution of purity in each detector is quite 
different from the evolution of entanglement between them. There is no "sudden death" of quantum entanglement 
between one detector and the rest, and the late-time purity is always less than one. This is a consequence of the 
direct interaction between each detector and the field. 



E. Disentanglement time and separation of detectors 

Naively one may think that spatially the closer Alice and Bob are, the larger ( TZa,TIb ), so the disentanglement 
time gets shorter. Our results in Fig. [8] show that this is not true. One can hardly see any simple relation between 
the initial separation of the two detectors and the disentanglement time. The larger the separation, the stronger the 
entanglement (S more negative) at some moments, but weaker at others. 

We have considered the case with both detectors being at rest and spatially separated. While the higher order 
corrections from mutual influences are complicated and quite hard to handle as t increases, for separation d large 
enough that the mutual influences cannot reach in time (namely, d > ct), we get a similar result that no simple 
proportionality exists between the separation and the degree of entanglement of the two detectors. As the system 
evolves, sometimes the entanglement is stronger for larger separation, sometimes it is weaker. 



F. How generic are the features contained in this model? 



Our results show that the disentanglement time in our model is finite in all cases considered. But how generic are 
these results? 

Our model has two UD detectors spatially separated in (3+l)D spacetime, one stationary and another uniformly 
accelerating and running away, without any direct interaction with each other. The presence of the event horizon for 
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FIG. 7: Purities of Alice (Va) and Rob (Vb)- (Top) Va and Vb in the ultraweak-coupling regime, with the same parameters 
as those in Fig. [2] In this regime Va and Vb each in its own proper time has virtually the same behavior. (Middle) Va and 
Vb in the strong-coupling regime, with the same parameters as those in Fig. [6] Again their difference is very tiny. (Bottom) 
Va and Vb in high acceleration limit, with the same parameters as those for the solid curve in Fig. [5] Now their difference is 
evident. 

Rob has the effect of curtailing the mutual influences higher than a certain order. When we focus on one single UD 
detector, it behaves like the QBM of a harmonic oscillator interacting with an Ohmic bath provided by the (3+l)D 
scalar quantum field [lj|. 

If there exists a direct interaction between the two detectors, there will certainly be residual entanglement between 
them at late times, just like the late-time entanglement between each detector and the field in our model. 

The scalar field in our model is living in a free space. The presence of boundaries will change the result. The 
backreaction started from each detector when reflected by boundaries will add another memory effect to the dynamics 
of entanglement and affect the late-time behavior of detectors. One example is an electromagnetic field in a perfect 
cavity. Spacetime dimension also matters. Different spacetime dimensions or boundaries give different spectral density 
functions of the field experienced by the detectors [lj|. Moreover, in (n + 1)D, n > 3 the spatial separation d implies 
a suppression of mutual influences by a factor l/d n ~ 2 , so the mutual influences would be negligible in the cases with 
large spatial separation. 

Spatial separation also contributes to the retardation of mutual influences, and since Alice is at rest and Rob is in 
uniform acceleration, one detector's dynamics is generally out of phase from the backreaction of the other. If the two 
detectors are located very close to each other through the whole dynamical process, or the two detectors are separated 
at a distance where resonance (d « 2irc/Q) is set up, the features of the disentanglement process can be different. For 
example, residual entanglement at late times under specific conditions in the 2HO QBM model has been reported in 
flU [lt| ; that model is equivalent to a detector theory with two UD detectors sitting at the same point in 3D space, 
when the HO bath is ohmic. 

Therefore, we believe that our results in this paper are generic for two well-separated detectors or atoms without 
direct interaction with each other, but coupled to a common quantum field in (n + 1)D, (n > 3) Minkowski space 
without boundary. 
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FIG. 8: (Upper) E with initial separation di ~ 0.5025 (b = 2.01) for the solid curve and di = 0.9999 (b = 10000) for the 
dashed curve. Alice is located at z\ = 1/6 and the initial separation between Alice and Rob at t = is di = a -1 — b^ 1 . Here 
a = 1, 7 = 0.025, Ao = Ai = 50, and other parameters are the same as those in Fig. [2] One can see that short-time revivals 
of entanglement occur during t « 10 ~ 11 and around t = 12. (Lower) The dependence of E on the initial separation is not 
evident systematically: from left to right are E at t = 0.5, 2.5, 3.2, respectively. 
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APPENDIX A: EXPRESSIONS FOR MODE FUNCTIONS 



For the cases with b > 2a, one has 
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q B + \r, k) = 6(t)^ jf dr'K (t T>)ti + \z B (T>), k) 

+ (r - n) ^ jf ^^n-i / ^'-^ " k). (A6) 

where n = a _1 ln(6/a), tj. = a" 1 - b~ x , t 2 = bo,- 2 - b^ 1 (see Fig. [J), r_(i') = a _1 lno (f + *-( r ') = 

— a _1 e~ ar , W± = l±[(Q r + i-f)/9l}, and if (x) = e~ JX sinflx. In ultraweak-coupling limit, we shall neglect O(Aq) 
terms. 



APPENDIX B: TWO DETECTORS AT REST WEAKLY COUPLED TO A THERMAL BATH 

In the ultraweak-coupling limit, for Q A and Q B both inertial and in contact with a thermal bath of temperature 
T, one has 



where j = A, B, and 
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Then the quantity E reads 



where 
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and Y n — I„| ^(sffl,Hs/«)' n = 0,2,4. For T > 0, if a 2 (3 2 h 2 , E < at t — and E > as t — > 00, so the 
disentanglement time tdE must be finite. It is easy to verify that there exists only one solution for E^s) = 0: For 
a 2 /? 2 > h 2 , 
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and the disentanglement time for a 2 /? 2 < ft 2 is tdE(X n — > Y n )\a?8 2 >fr 2 - At high temperature limit T — > oo, we have 

tdE ~ 
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